Langevin description of speckle dynamics 
in nonlinear disordered media 
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We formulate a Langevin description of dynamics of a speckle pattern resulting from the multiple 
scattering of a coherent wave in a nonlinear disordered medium. The speckle pattern exhibits 
instability with respect to periodic excitations at frequencies Q below some f2 ma x, provided that 
the nonlinearity exceeds some ^-dependent threshold. A transition of the speckle pattern from a 
stationary state to the chaotic evolution is predicted upon increasing nonlinearity. The shortest 
typical time scale of chaotic intensity fluctuations is of the order of l/fimax- 



Propagation of a coherent wave in a disordered medium 
is diffusive if A -C I -C L, where A is the wavelength, I 
is the mean free path, and L is the size of the medium 
JD. While the wave undergoes multiple scattering and 
the spatial distribution of the scattered intensity looks 
quite irregular (speckle pattern), the coherence of the 
wave is not destroyed and various coherent phenomena 
can be observed: enhanced backscattering, short- and 
long-range intensity correlations, universal conductance 
fluctuations, etc. (see Rcfs. [|, || for reviews). Available 
studies of nonlinear phenomena for diffuse waves include 
calculations of the enhanced backscattering cone at fun- 
damental Q and doubled || frequencies, investigations 
of optical phase conjugation ||, studies of correlations 
in transmission and reflection coefficients of the second 
harmonic and fundamental Q waves, an extension of 
the standard diagrammatic technique to nonlinear disor- 
dered media ||] , and a study of persistent hole burning 
in multiple-scattering media |To|| . 

After realizing that the sensitivity of the speckle pat- 
tern to changes of the scattering potential diverges for a 
sufficiently strong nonlinearity |11| , a new phenomenon, 
the temporal instability of the multiple-scattering speckle 
pattern in a disordered medium with cubic nonlinearity, 
has recently been predicted The speckle pattern 

is expected to become unstable and to exhibit sponta- 
neous fluctuations if the nonlinearity exceeds some crit- 
ical value. Although of primary importance in view of 
the possible experimental observation of the instability 
phenomenon, the dynamics of spontaneous intensity fluc- 
tuations, their nature and associated characteristic time 
scales have not yet been studied up to now. 

In the present paper we formulate the dynamic 
Langevin description of spontaneous intensity fluctua- 
tions in a nonlinear disordered medium. Our theoretical 
method can be viewed as an extension of the station- 
ary Langevin approach introduced in Ref. jOJ the latter 
being inadequate to describe the dynamics of speckles. 
Analysis of the speckle pattern stability with respect to 



weak periodic excitations shows that if the effective non- 
linearity parameter p — An 2 (L/l) 3 exceeds some criti- 
cal value p c ~ I (where An is the typical value of the 
nonlinear correction to the refractive index) the speckle 
pattern becomes unstable with respect to periodic exci- 
tations at frequencies inside some limited low-frequency 
interval, and the maximal Lyapunov exponent becomes 
positive. This allows us to describe the chaotic nature 
of spontaneous intensity fluctuations beyond the abso- 
lute instability threshold p = p c and to estimate their 
characteristic time scale. 

We consider a scalar wave propagating in a nonlinear 
disordered medium and described by the following wave 
equation: 



2 r 



c 2 dt 2 
= J{r,t), 



l + e(r,<)+e 2 |^(r,t)| 2 U(r,t) 



(1) 
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where J(r,t) = Jo(r) exp(— iwot) is a monochromatic 
source term, c denotes the speed of wave in the aver- 
age medium, e(r,t) is the fractional fluctuation of the 
dielectric constant at frequency uj$ (possibly slowly vary- 
ing in time), and £2 is a nonlinear constant. Equation 
(|l|) describes, e.g., propagation of optical waves in me- 
dia with intensity-dependent refractive index Jl^| in the 
scalar approximation and neglecting the generation of the 
third optical harmonics. The latter assumption is justi- 
fied in the absence of phase matching |ll| or, more pre- 
cisely, when |/c(3wo) — 3k(wo)\ I ^ 1> where k(ui) is the 
wavenumber at frequency uj. 

Consider first a linear medium (e 2 = 0) of typical size 
L and a white-noise Gaussian disorder: (e(r,t)e(ri,t)) = 
Air/ (fcoO^( r— r i)j where k = k(uj ) = uj /c. Let the time 
variations of e(r, t) be random, stationary, and arbitrary 
slow, so that the time scale of the resulting variations 
of the amplitude <p(r,t) of i>(r,t) — tp(r,t) exp(—iwot) 
is much larger than the typical time between two suc- 
cessive scattering events l/c. For L 3> I and far enough 
from the boundaries of the disordered sample, the av- 
erage intensity (I(r)) then obeys the diffusion equation 
JTij , while the long-range correlation of intensity fluctu- 
ations SI(r,t) = I(r,t) — (I(r)) can be found by solving 
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FIG. 1: Diagrams contributing to Eqs. and (^). Solid 
lines denote the wave field ip and the complex conjugated 
field ip* . Dashed lines denote scattering of ip and i/i* on the 
same heterogeneity. The diagrams (b) and (c) are obtained 
by inserting two fcg vertices (denoted by wavy lines) to the 
diagram (a) at (r',t') and (r{,t{), respectively. 
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functions q(r, r ', At = t — t') = <5j oxt (r, t)/Se(r ' t ') can 
be found by a functional differentiation of Eq. (g): 




the Langevin equation J15[ : 



d 



j- t 5I(r, t) - DV 2 5I(v, t) = -V • j ext (r, t), 



(2) 



where I{r,t) — \ip(r, i)\ , £> = cZ/3 is the diffusion con- 
stant, and jext(r, t) are random external Langevin cur- 
rents: 

i£kr,t)£kri,*i) 
= 27rZc 2 /(3fc 2 ) I (^(r, i)^ (r, ij)) | 2 %«5(r - n). (3) 

The diagram corresponding to Eq. (||) is shown in Fig 
0(a). 

For a given Jo(r), the current j ext (r,t) is a "finger- 
print" of the disorder e(r, t). An infinitesimal variation 
Ae(r, t) of the dielectric constant will modify j cx t(r, t) by 
a small amount 

Ajext(r,i) 

= f d 3 r' / di'q(r,r',*-*') Ae( r '»*')> ( 4 ) 

where the spatial integral is over the volume V of the 
sample, we neglect the terms of the second and higher 
orders in Ae(r, t), and the correlation of random response 



(r.r'.AtJgWCri.ri'.Atx) 



= 3^D 2 (c 2 /Z) ( 5. i ^(r-r 1 ) 

x [<I(r')> G(r',r{; At - A^G^.r; At x ) <J(r)> 
+ (7( ri ')) G(r x ', r'; Afc - Ai)G(r', r; At) (7(r)) 
-(J(r')>G(r',r;At)(J(r 1 '))GW,r;At 1 )]. (5) 

Here |r — r'|, |r — r{|, |r' — r{| ^ Z is assumed and 
G(r, ri; At) is the Green's function of Eq. (|J). The dia- 
grams contributing to Eq. (j^) are shown in Fig. [l](b, c). 
In the stationary limit Ae(r, t) = Ae(r), Eqs. (||) and (||) 
reduce to Eqs. (6) and (7) of Ref. 

The Langevin description of intensity fluctuations in 
disordered media can be extended to the nonlinear case 
(e2 7^ 0). To this end, we consider time-independent e: 
e(r,t) = e(r), and assume that the diffusion constant D 
and the mean free path Z are not affected by the nonlin- 
earity. The latter assumption is valid if An 2 kol <C 1 [ pd| , 
where An = (s 2 /2)io 1S the typical value of the nonlin- 
ear correction to the refractive index and Iq ~ {I( r )) is 
the typical value of the average intensity in the medium. 
We now admit that in a nonlinear medium, the total di- 
electric constant contains a linear contribution 1 + e(r) 
that we assumed to be time-independent, and a nonlinear 
contribution £2^(r, i) that can vary with time. The varia- 
tion of the total dielectric constant can be therefore only 
due its nonlinear part, and we can identify the infinitesi- 
mal variation of the dielectric constant Ae(r',t') in Eq. 
(g) with e 2 AI(r',t / ), where AJ(r',t') is the change of 
the intensity at r ' during some infinitesimal time interval 
(t',t' + 5t): AI(r',t')=I(r',t' + 5t)-I(r',t'). Substi- 
tuting Ae(r',i') = e 2 AI(r',t') into Eq. (g), noting that 
Ajext(r, t) = j ex t(r, t + St)- jext(r, t), and dividing both 
sides of the resulting equation by 5t —> 0, we obtain the 
following dynamic equation: 



d_ 

dt 



j ext (r,t) = £ 2 / d 3 r' f dAtq(r,r',At) 
Jv Jo 



x -5I{v',t-At), 



(6) 



where we make use of the fact that I(r, t) = (I(r)) + 
51{r,t) and hence AI(r,t) = A[5I(r,t)}, since (I(r)) is 
time-independent . 

Equations (||) and (^) form a self-consistent system of 
equations: Eq. (^) governs the spatio-temporal evolution 
of the intensity fluctuations <5/(r, t) due to the Langevin 
currents j ex t(r, t), while Eq. (||) describes the distributed 
feedback mechanism, leading to variations of j cxt (r,t) de- 
pending on the changes of 51 (r, t). Note that Eq. (|J) is a 
linearized equation: only the terms linear in the nonlinear 
contribution to the dielectric constant £2-T(r, t) ar e kept, 
which is justified as long as s^Iir, t) <C 1. In certain cir- 
cumstances (see below) , the linearized nature of Eq. (||) 
may result in the exponential growth of its solution with 
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FIG. 2: Surface describing the stability of the multiple- 
scattering speckle pattern in a nonlinear disordered medium 
with open boundaries. At given effective nonlinearity param- 
eter p and frequency Q, the surface defines the Lyapunov ex- 
ponent A. If A > 0, the speckle pattern is unstable with 
respect to periodic excitations at frequency Q. 

time, and in this sense Eqs. (||) and (JsJ) are analogous 
to the equations of linear stability analysis commonly 
used to study the stability of nonlinear systems (see, e.g., 
Rcfs. [l6| and |l7| for examples of nonlinear optical systems 
exhibiting instabilities). Hence, although Eqs. (||) and 
@ allow us to study the stability of the speckle pattern 
and the characteristic time scales of spontaneous inten- 
sity fluctuations beyond the instability threshold, they 
cannot be used to determine the amplitude of these fluc- 
tuations. 

Consider now an infinitesimal periodic excitation of 
the static speckle pattern: 6I(r,t) — SI(r, v) exp(iisb), 
where v = il — iA ^ and fl > 0. Such an excita- 
tion can be either damped or amplified, depending on 
the sign of the Lyapunov exponent A. The value of A 
is determined by two competing processes: on the one 
hand, diffusion tends to smear the excitation out, while 
on the other hand, the distributed feedback sustains its 
existence. The mathematical description of this compe- 
tition is provided by Eqs. (|^) and (^), that after the sub- 
stitution of 6I(r,t) = SI(r, v) exp(ivt) [and similarly for 
j ex t(r, t)] lead to the following equation (see Appendix |A] 
for the details of derivations): 

p~F(Q/fl D ,A/fl D ). (7) 

Here p — An 2 (L/l) 3 is the effective nonlinearity param- 
eter, the function F is shown in Fig. ^|, and a numeri- 
cal factor of order unity is omitted. To obtain Eq. (Q), 
we have assumed the disordered sample to have open 
boundaries (i.e. the diffusing wave leaves the sample 
when it reaches a boundary) and have taken the limits 
of large sample size (L/l 3> kol) and moderate frequency 
fl <C flD[L/(k l 2 )] 2 , where Qd = D/L 2 is the inverse of 
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FIG. 3: Main plot: frequency-dependent "phase diagram" 
of the multiple-scattering speckle pattern in a nonlinear dis- 
ordered medium with open (solid line) or reflecting (dashed 
line) boundaries. For a given f2, p should exceed the plotted 
threshold value for the instability to develop. Inset: maximal 
Lyapunov exponent as a function of the effective nonlinearity 
parameter p. The dotted lines show A max = and p = 1. 



the typical time needed for a multiple-scattered wave to 
diffuse through the disordered sample. 

It follows from Fig. || that for a given frequency SI the 
sign of the Lyapunov exponent A depends on the value of 
p. Excitations at frequencies SI corresponding to A < 
are damped exponentially and thus soon disappear. In 
contrast, excitations at frequencies S7 corresponding to 
A > are exponentially amplified, which signifies the in- 
stability of the speckle pattern with respect to excitations 
at such frequencies. Noting that A is always negative for 
P < 1, we conclude that all excitation are damped in 
this case and the speckle pattern is absolutely stable. In 
an experiment, any spontaneous excitation of the static 
speckle pattern will be suppressed and the speckle pat- 
tern will be independent of time: SI(r,t) = 51 (r), as 
in the linear case. When p > I, an interval of frequen- 
cies < SI < Sl max starts to open up with A > 0. The 
speckle pattern thus becomes unstable with respect to 
excitations at low frequencies. In an experiment, any 
spontaneous excitation of the static speckle pattern at 
frequency SI G (0, S7 max ) will be amplified and one will 
observe a time- varying speckle pattern SI(r,t). 

The border between stable (A < 0) and unstable 
(A > 0) regimes is shown in Fig. |3| by a solid line. The 
instability threshold increases with SI. At SI <C the 
exact functional dependence of the threshold on SI is 
rather sensitive to the peculiarities of the disordered sam- 
ple (e.g., its geometry and conditions on the boundaries), 
since for such slow oscillations the feedback mechanism 
is ensured by partial waves that have long path lengths 
s > L 2 1 1 and hence "feel" the presence of the bound- 
aries and the shape of the sample. Using the analytic 
expression for the function F (Q/Qd, A/ilp) derived in 
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Appendix we find p - 1 - (ft/^) 1 /" with a ~ 1/2. 
This yields f2 max ~ ^d(p — 1)", and the shortest typi- 
cal time scale r of spontaneous intensity fluctuations can 
be estimated as r ~ f2~ ax ~ ^^{p — 1)~ Q . At high fre- 
quencies Q £Id we find p~ (fi/fi^) 1 / 2 , Sl max ~ £IdP 2 , 
and r ~ p~ 2 , respectively. The latter results, on the 
contrary, are weakly sensitive to the peculiarities of the 
sample, since for the fast oscillations the feedback mech- 
anism is due to relatively short diffusion paths that do 
not reach the boundaries of the sample. 

The rise of the instability threshold with f2 can be qual- 
itatively understood by considering the phase difference 
A0(At) between two waves traveling through the disor- 
dered sample along the same diffusion path but separated 
in time by At ~ f2 -1 . If e(r, t) changes slowly with time, 
A<fi(At) comprises two contributions: A4>L,{At), which is 
the phase difference in a linear medium, and A</>jvx(A£), 
which is the additional phase difference due to the non- 
linearity. The second moment of the latter is 



<A^ L (At)) 



^o e 2 / dsi / ds 2 
(A/(ri,At)AJ(r a ,At)>, 



(8) 



where Sj is a curvilinear coordinate of the point rj, 
the integrals are along the diffusion path of typ- 
ical length s ~ L 2 /l, and AJ(ri,At) denotes the 
change of the intensity at during the time At. For 
0, <c Q,d we can assume that (Ai Yri , At)AJ(r2, At)) ~ 
(A0|(At)> (*/(n,0)W(r 2 ,0)) @. Taking 
(A(j) 2 NL (At)) > (A<f>l(At)) to be the instability 
condition for the excitation of the speckle pattern at 
frequency Q, and noting that (<$/(ri, 0)<5/(r 2 , 0)} ~ 
IqKUqI \y\ — r 2 |) ]l5fl , we recover p > 1 as the instability 
criterion. If Q 3> the long-range intensity correla- 
tion establishes only for ri — r 2 < (DAt) 1 / 2 , and the 
instability condition becomes p > (fi/f^) 1 / 2 ^> 1. 

The positive sign of the maximal Lyapunov exponent 
A max for p > 1 (solid line in the inset of Fig. [|), as well 
as the continuous spectrum < < O max of frequencies 
with A > 0, are hallmarks of chaotic behavior |Q. A 
sharp transition of the speckle pattern to chaos at p = 1 is 
reminiscent of the behavior observed in nonlinear systems 
with large (infinite) number of degrees of freedom (e.g., 
random neural networks with an infinitely large number 
of nodes p9|) and should be contrasted to the "route to 
chaos" through a sequence of bifurcations, characteristic 
of low-dimensional nonlinear systems [ jig ]. As one can 
see from Fig. |^, the scaling of A max with p — 1 appears 
to be roughly linear for \p — 1| <C 1: A max ~ VLd{p — l)* 3 
with (3 ~ 1. 

To demonstrate the sensitivity of the results obtained 
for <C £Id to the peculiarities of the disordered sam- 
ple, we briefly consider a sample with reflecting bound- 
aries (dashed lines in Fig. ||). All calculations can be 
carried out in the same way as for the sample with open 
boundaries (see Appendix 0) , assuming that the Green's 
function of Eq. (0) is approximately the same as in the 



infinite medium: G(r,ri; At) ~ Gq(t, ri; At). We find 
that the absolute instability threshold p c is roughly 2 
times lower than in the open geometry, a ~ (3 ~ 2, and 
A max = for p < p c . For an arbitrary sample of disor- 
dered nonlinear medium we expect O max ~ £Id{p — Pc) a 
and A max ~ £Id{p - Pc) 13 for p - p c < 1 and p > p c , 
where p c ~ 1, 1/2 < a < 2, and 1 < (3 < 2. By analogy 
p8| with the theory of phase transitions, A max and /3 can 
be identified with the order parameter and the critical 
exponent, respectively. 

Finally, it is worthwhile to note that Eqs. (|J) and 
(^) can also be derived from a time-dependent disor- 
dered nonlinear Schrodinger equation with a potential 
u(r) + g \ip(r, t)\ 2 . Upon the substitutions uj — > E/h, 
kl -> 2mE/h 2 , e(r) -> [-tt(r)/£], and e 2 -> (-.g/S) 
[where E is the energy of the incident Schrodinger wave, 
m is the particle mass, u(r) is the disordered potential, 
and g is the nonlinear constant], our analysis is there- 
fore valid in this case too. The analogy between the 
wave equation (Q) and the Schrodinger equation is known 
for the stationary case, when the solution ip(r,t) can be 
represented as ipo(r) exp(— iuj^t). However, the dynamic 
solutions of the two equations differ due to different dis- 
persion relations. Although the present paper deals with 
dynamic speckle patterns, their temporal fluctuations are 
assumed to be slow and the analogy between the wave 
and Schrodinger equations is recovered within the accu- 
racy of our analysis. 
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APPENDIX A: DERIVATION OF EQ. (g) 

In this Appendix we provide a derivation of Eq. 
(0) from Eqs. (§) and (§). Substituting SI(r,t) = 
<J/(r, z/) exp(iz4) and j ext (r,i) = jext(r, v) exp(ii4) into 
the two latter equations, we obtain 

i V 5I(v,y)-DV 2 SI{v, V ) = -V-jextM, (Al) 

^jext(r,z/) =ive 2 / dV / dAtq(r,r', At) 
Jv Jo 

(A2) 



xSI(r', v) exp(— ivAt) 



If v = 0, Eq. (A2) is trivial and the statistical proper- 
ties of j ex t(r,0) are determined by Eq. (|J) with t = ti, 
the same equation as in the case of the linear medium, 
while the static, time-independent part of the intensity 
fluctuation 6I(r, 0) is found by solving the stationary 
Langevin equation [Eq. QA1| ) with v = 0]. Hence, the 
time- independent part of the speckle pattern remains the 
same as in the linear medium. If, in contrast, v ^ (as 
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we a ssume in the main text), we divide both sides of Eq. 
( |A2| ) by iv, multiply the i-th Cartesian component of the 
resulting equation by the j-th Cartesian component of 
a similar equation for j* xt (ri,i/), and average the result 
over disorder. This yields 



j«(r,i/)j£> (rijI/) 



d 3 r' 



d 3 r{ (SI(r',u)SI*(r[,u)) 



poo POO 

x dAt dAii (gWfr.r'.AtJgW'Crx.r^Ati 
Jo Jo ' 

x exp[-ivAt + iv*Ah)], (A3) 



where j^ t (r,i/) denotes the i-th Cartesian component 
of jext(r, v). After the substitution of Eq. (g) for 
(q^(r,r', At)q^*(ri, r{, Ati)), the time integrations in 
Eq. © yield 

rOO POO 

J dAtjf dAh (^(r^At^fa.r^Ati: 

x exp[— wAt + w*Ati)] 
= 3 7 rL> 2 (c 2 /0' 5 'u^( I '- r i) 



x[(7(r'))G(r' 1 r 1 V)G*(rl,r; i ,-^)(7(r)) 
+ (J(rO)G*(r 1 / 1 r';i/)G(r / J r;i/-i/*)(J(r)> 
-<I(r')}G(r',r;i/) (1(4)) G* (r[, r; i/)] , 



(A4) 



where G(r, r^; !/) is the Fourier transform of G(r, ri; At). 
Equation ( [A3] ) can now be rewritten as 

j&(r, i/)) = A(r,i/)^5(r- ri ),(A5) 



where 



37rU 2 (c 2 /0£2 / d 3 r' / d 3 r{ 
Jv Jv 



x(8I(v',v)8I*(v[,v)) 

x[(I(r'))G(r',r{;v)G*(r{,r;v-v*)(I(r)) 
+ (I(v[))G*(v' L ,v';v)G(v',v;v-v*)(I(v)) 
-(/(r'))G(r',r;^)(/(rO)G*W,r;^)]. (A6) 

In the following, we replace both (J(r)) and A(r, v) by 
their spatial averages 1$ and A(z/), respectively. This sim- 
plifies the further analysis considerably, while can only af- 
fect the final result by a numerical factor of order unity, 
since (1(f)) and A(r, v) do not change significantly as 
long as the point r is far enough from the sample bound- 



We now admit that Eq. (A5) for the correlation func- 
tion of Langevin currents at v ^ in a nonlinear 
medium has a form similar to Eq. (|J) for Langevin 
currents in a linear mediu m. This allows us to pro- 
ceed with analysis of Eq. ( |Al[ ) in the same way as it 
was done for Eq. (Q) in the linear medium |n|. To 
simplify further calculations, we assume that the dis- 
ordered sample has open boundaries (i.e. that the 



multiple-scattered waves leave the sample when they 
reach the boundary) and hence the Green's function 
of Eq. (||), G(r, ri;Ai), can be approximately writ- 
ten as Go(r, ri; At) exp(— Sip At), where Go(r,ri; At) = 
(iirDAt)-^ 2 exp[- |r - r x \ 2 /(AD At)] is the Green's 
function in the infinite medium, exp(— flu At) describes 
the leakage of the wave through the sample bounda ries , 
and H.D = D/L 2 . We now write the solution of Eq. ( Al) 
as 

5I(v,v) = - [ d 3 r'G(r,r»[V-j oxt (r»] 
Jv 

= / d 3 r'[VG(r,r»-j oxt (r»], (A7) 
Jv 

where the second line is obtained as a result of integra- 
tion by parts, assuming G(r, r '; v)j^ t (r', v) = at the 
boundary of the disordered sample. Multiplying Eq. (A7) 
by a similar equation for 8I*(tujA, performing the aver- 
aging over disorder using Eq. (A5), and carrying out the 
necessary integrations, we obtain 

(5I(T,v)5r( ri , V )) 



My) 

D 



ReG(r, n; v) 



Im^ 

Rcj/ 



ImG(r, T\\v) 



.(AS) 



Substituting Eq. (|A|) into Eq. QAq), dividing both 
sides of the resulting equation by A(u) ^ 0, recalling 
that v = Q — iA, and performing a change of variables 
R = [n/(2D)] 1 / 2 v (and similarly for r' and r{), we ob- 
tain 



1 = CiAn 2 (L/lfh(Il,n/VL D ,AID, D ), 



(A9) 



where C\ is a numerical constant, Art = (e 2 /2)-Zo; an d 
the dimensionless function h is defined as 

h(K,n/n D ,A/n D ) = (n D /n) 3/2 J d 3 R' J d z R[ 

x[g(R',R{; 1 )G* 1 (R[,R;ii) 
+ g*(R[,R'; 1 )g 1 (R',R; ll ) 
-g(R',R; 1 )g*(R[,R; 1 )} 



Reg(R',R 1 '; 7 ) + ^Im0(R', R{; 7 ) 



(A10) 



C?(R,Ri;7) 
1 

~ |R-Ri| 

5i (R, Ri;7i 
1 



exp[-( 7 + V7)|R-Ri|], (AH) 



exp(-2 7 i |R-Ri| 

1/2 




(A12) 



(A13) 
(A14) 
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We now assume that the disordered sample has 
the shape of a sphere centered at the origin 
and that h(0, fi/firj, A/Od) provides a good esti- 
mation of h for the points R located far enough 
from the boundaries. Defining F(f2/f2.D, A/Qd) = 
h(0, O,0)//i(0, Q/Qd, A/Qd) and introducing the effec- 
tive non linearity parameter p = An 2 (L/l) 3 , we rewrite 
Eq. {M) as 



p = c 2 F(si/n D , A/n D ), 



(A15) 



where C2 is a numerical factor of order unity. Since we 
have already made some approximations that affect the 
final result by a numerical factor of order u nity (e.g., 
we replaced (/(r)) by I ), we omit C2 in Eq. ( A15 ) and 
obtain Eq. (fiT) of the main text. The most of integrations 



in Eq. ( Aldf ) can be performed analytically, while the 
remaining integrations are easily carried out numerically, 



allowing us to de term ine the value of A for given p, fir), 
and Q from Eq. ( Af5 ). 



A n imp ortant comment is in order in connection with 
Eq. ( Af 5 ) and the analysis it results from. The correla- 
tion function of intensity fluctu atio ns (SI(r, v)SI* (ri, v)) 
entering into Eqs. (A3) and (A6) contains, in princi- 
ple, not only the long-range contribution given by Eq. 
(|A§|), but also a short-range one {5I(r, u)5I*(ri, v)) ~ 
(l/k^I^S^ — ri). The latter contribution has been 
neglected in our analysis, which is justified for large 
enough sample size (L/l 3> fcoO and moderate frequency 
il <C fin[£/(fco/ 2 )] 2 . If one of the above inequalities is 
violated, the roles played by the short- and long-range 
contributions to the correlation function of intensity fluc- 
tuations in development of the instability become com- 
parable, and the above analysis is no longer valid. 
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